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ABSTRACT 


A theoretical investigation of gas flow inside a multilayer insulation system 
has been made for the case of the broadside pumping process. A set of simultaneous 
firsh-order differential equations for the temperature and pressure of the gas mixture 
was obtained by considering the diffusion mechanism of the gas molecules through 
the perforations on the insulation layers. A modified Runge-Kutta method was used 
for numerical experiment. The numerical stability problem was investigated. It 
has been shown that when the relaxation time is small compared with the time period 
over which the gas properties change appreciably, the set of differential equations 
can be replaced by a set of algebraic equations for solution. Numerical examples 
were given and comparisons with experimental data were made. 
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Subscri pts: 


Superscripts: 


area of one insulation layer (one side) 

total area of perforation on one insulation layer 

a 1 - (1 - c) A 

perforation coefficient c = a/A 
Boltzmann constant 

number of interstitial spaces between insulation layers 
number of kinds of absorbed gases plus 1 
mass of a gas molecule or an integer 
number of gas molecules 

number density of gas molecules or an integer 
pressure 

outgassing rates, see Equations (4), (7), and (12) 

temperature 

time 

volume 

mean thermal velocity of gas molecules 
number of degrees of freedom of a gas molecule 
relaxation time 

subscript c refers to quantities inside the vacuum chamber, f 
to condition ot the fuel tank, a to atmospheric condition, othe 
subscripts refer to quantities in a certain interstitial space 
superscript 1 refers to the purge gas,' other superscripts refer 
to a certain kind of gas molecules 



1. Introduction 


Recent developments in propulsion technology have stimulated interest in 
the study of multi layer insulation systems for the fuel tank of a rocket booster. The 
purpose of multilayer insulations is to eliminate heat conduction between the fuel 
tank and the environment. It is then desirable that the spaces between the insula- 
tion layers be maintained at low pressures. However, it is not practical to require 
complete vacuum in the spaces between the insulation layers, although it would be 
ideal in principle, for the reasons that complete vacuum is hardly attainable and 

that when the pressure gets too low, deformation of the layers would cause solid 

1 

contacts, which in turn would induce heat conduction. It is therefore customary 

-4 

in practice to require a pressure to be of the order of 10 torr. This low pressure 

is usually achieved before launch from the ambient pressure by some pumping device. 

There have been two pumping procedures in practice; namely, the broadside pumping 

by which the direction of gas flow is perpendicular to the insulation layers, and 

the edge pumping which causesthe gas to flow parallel to the layers. Experiments 
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of both evacuation procedures have been performed and results reported. ' The 
present work is a theoretical analysis for the broadside pumping process. 

In the case of broadside pumping, the insulation layers are perforated so 
that gas molecules can go through small holes on the layers, resulting in streaming 
gas motion perpendicular to the insulation layers. It is obvious that, in order to 
minimize layer deformation and solid contacts, small pressure differentials across 
the layers should be maintained during evacuation. This can be achieved by small 
perforation, i.e., the total area of the holes on an insulation layer is small com- 
pared to the total area of that layer. If we further require that the diameter of a 
single holebesmall or comparable to the distance between two successive layers, 
there is the advantage of uniformity of gas motion, which renders convenience dur- 
ing operation. 

** A 

A multilayer insulation system consists of a large number On the order of 10; of 
extremely thin sheets of low thermal conductivity and lightweight materials, kept 
parallel to one another with a total thickness of about 1". 
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2. Formulation 


Consider a fuel tank of simple geometry wrapped with multilayer insulation 
and situated inside a vacuum chamber <> Consider the broadside pumping process 
and assume that the perforation areas are small . Then, a gas molecule in any inter- 
stitial space between two successive layers will, on the average, collide many times 
with the walls or with some other gas molecules in that interstice before getting 
through a hole to a neighboring interstice. Thus, it is plaussible to assume that the 
gas in any interstitial space is in thermal equilibrium with the temperature and 
pressure pertinent to that interstice, and the motion of the gas is simply a diffusion 
process . 

Consider the gas inside the insulation system as a mixture of a certain kind 
of purge gas and a number of different kinds of gases originally absorbed in the in- 
sulation materials. Let be the total area of the outermost insulation layer (one 
side), and c^ the perforation coefficient of the same layer (the ratio of the total 
perforation area on the first layer to A^), and so on. Let be the volume of 
the interstice between the first and the second insulation layers, and the number 
of molecules in V^, V the volume of the interstice between the second and the 
third insulation layers, and ^ the number of molecules in V^, and so on. Since 
the gases in any interstice are in thermal equilibrium, each component gas has a 
Maxwellian distribution with its local density and the common local temperature. 

From the kinetic theory of gases, we have, for the gases in the i th interstice. 
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where a. - c. A. , a 1 . ” (1 “ c.) A. , k is the Boltzmann constant, 6 the num- 
ber of degrees of freedom of motion of a gas molecule, n., v. and T. are the number 

i i i 

density, the mean thermal velocity and the temperature of the gas molecules in the 
i th interstice, respectively. The superscript j refers to the purge gas if j = 1, 
or the absorbed gas of kind j otherwise. The number of interstitial spaces is L, 
and the number of kinds of the absorbed gases is (M-l) so that 



Finally, 



0 (j = 1) 


(4) 


| number of outgassing molecules, of kind j from the 
walls binding the i th interstice, per unit area and 
time ( j ^ 1 ). 

Equations (1) state the conservation of mass and Eqs. (2) the conservation of energy. 
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Since the L th Interstice is nect to the fuel tank, we may assume that 

a = 0 and T = T , where T, is the temperature of the fuel tank. Since inter- 
L + 1 L f f 

stice 1 is the outermost, denoting the number densities, the mean thermal velocities 

/• \ /• \ 

and the temperature of the molecules in the vacuum chamber by n^ , v an< ^ 

T , respectively, we have n ^ = n v ^ = v ^ and T ~T • Notice that 
as the volume of the vacuum chamber is very large, n (j ^ 1) are vanishingly 

/ 1 \ C /• \ 

small compared to n ^ "at ail times and we may assume that n ^ = O(j^l). 

C (i ) c 

By making this assumption, v (j ^ 1 ) become meaningless and they drop out from 

o 

Eqs. (l)and (2) automatically . Therefore, in Eqs. (l)and (2), we use 


a. , . = 0 , T. = T_, T =T , v 0) = v (1) 
L+l L foco c 


( 5 ) 


ondn G>. 

o c o 


0 (j^l) • 


Following p. ~n. k T. , where p. is the pressure in the ith interstice. 


we may assume that 


n.^ k T, so that 


i i i 
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which is plaussible when the pressure is not too high. 

Using Eqs. (3) and (6), and if the outgassing rates are expressed in the 
customary units of pressure times volume per unit area and time denoted by 

Q.® = kT. q.G* (7) 

i i ~i ' ' 

then, Eqs. (1), (2), and (5) can be written as 
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where p ^ is the chamber pressure due to the purge gas. 

G 

The mean thermal velocities are given by 




1/2 

,1=0, 1, , L; j = 1, 2, ,M . 
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where m y is the mass of a molecule of the j th kind. 

To obtain a solution to the problem, the oufgassing rates Q. ^ must be 

1 4 7 ^ 

known a priori, usually determined experimentally. ' Outgassing is a process 

in which gas molecules originally adsorbed or absorbed in a solid material leave 

that material under reduced pressure or elevated temperature. The outgassing rate 

of the insulation layers depends on the material and preconditioning of the layers, 

*** 

the temperature, pressure, time and the kinds of absorbed gases. For a specific 
material and preconditioning and a specific absorbed gas, Q.^ = Q.^ (f, p, t). 
Around the room temperature range, there is not outgassing at atmospheric or higher 
pressures. Outgassing occurs when the background pressure is reduced considerably 
below the atmospheric pressure. The outgassing rate increases as the pressure de- 
creases, and attains apprecible values only at very low pressures. However, in 
very low. pressure ranges, the outgassing rate is a weak function of pressure. There- 
fore, its dependence on the pressure may be approximated by 




G) 


( T, 0 


02 ) 


M 

= £ 
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where p. p. VJ/ , p^ is the atmospheric pressure and R^ are the out- 

gassing rates at extremely low pressures. 


*** 


The absorbed gases are mainly water vapor, with small amounts of and 
N 2 (Ref. 7). 
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The initial conditions associated with the system of Equations (8) and (9) 
are prescribed according to a specific situation encountered. Usually the evacuation 
process starts with atmospheric pressure inside the insulation system. In this case, 
the initial conditions can be written as 

P. (,) (o) - P a , P. (i) (o)= 0 (j / 1 ), i = 1, 2, 

T. (0) = T (i), i = 1, 2, , L - 1. 

Now, with the (M + 1) L - 1 initial conditions given by Eqs. (13) or pre- 
scribed otherwise and the outgassing rates Q. ^ known, and using Eqs. (10) and 
(11), Eqs. (8) and (9) constitute an appropriate set of (M + 1) L - 1 simultaneous 
first order differential equations for the same number of unknown functions of time, 
namely, T. (t), i = 1, 2, L - 1 and p. ^ (t), i = 1, , L; j = 1, ,M. 

In some engineering applications, e.g., the Apollo Telescope Mount insula- 
tion, the temperature is quite uniform across the entire insulation system, and the 
gas flow inside the system is nearly isothermal . In this case, the problem reduces 
to the solution of Eqs. (8) with Eqs. (10) and (1 1), and the initial conditions (13), 
when all temperatures are set equal to a constant. 


L 

(13) 


1 



3. Numerical Analysis 

In a multilayer insulation system, the number L usually is of the order of 

2 

10 . Thus, the governing Eqs. (8) and (9) consist of a large number of first order 
differential equations which are non-linear and coupled. The analytical solution 
of this system of equations can hardly be obtained, one therefore uses a direct nu- 
merical method for its solution. 

For simplicity in numerical experiments, we took an isothermal case, and 

assumed water vapor was the only absorbed gas present ( M = 2 ) . Furthermore, 

( 2 ) 

a hypothetical R v (T , t)and chamber conditions were used. When the temperature 

is constant, the energy equations (Eqs. (9) )are not needed, and the problem reduces 

to the solution of Eqs. (8) with initial conditions (13). In the process of numerical 

computation, however, the problem of instability occurred. A numerical program 

8 

using a modified Runge-Kutta method has been set up for this initial value problem. 
Computer experiments on this scheme showed that the system was stable when the 


values of T. 


0 ) = 


(4 V. /a. ) (it m ^/8 k T ) , which have a dimension of 


1/2 


the time, were iarge. However, in order to establish numerical stability for small 

values of Tj , the step sizes of integration had to be so small that a solution could 

not be obtained with a reasonable length of computer time. Since in ordinary ap- 

/• \ 

plications the values of T ^ are very small, we need to seek some alternative or 
approximation to the system of Eqs. (8), (9) and (13). 

For the case j = 1, Eqs. (8) can be written as 
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(14) 


where t. - t (a., T. ) - ( 4 V. / a. ) ( TTm ^ ^ / 8 k T. ) ^ 2 , t. _ i - t(q., T._j) 
T | "H T ^°i + V ^i + 1) and T. = T (a. + j, T.) have a dimension of the time. 


i 


and where r. _ ^ ~ T. / T._^ and r. + j - T. / T. + ^ are dimensionless quantities 
or order 1 . 
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Let the larqest of all the four t's defined above be T . Define the non- 

m 

dimensional temperatures and pressures as T'. = T. /T r , where T is a reference 

1 1 0) r 

temperature (say the room temperature), and p 1 . = p // p^. Let be the time 

period over which T\ and p.' change appreciably, and define the non-dimensional 

time as t* = t/T so that d T'. / dt 1 and dp 1 . / dt ' have their magnitudes of order 
e i r i 

1. 

In terms of T 1 ., p 1 . and t 1 , Eqs. (14) can be written as 
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dp r 

dt' 


P - 


d T'. 
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d t‘ 


T 'i [ r i-i s i-l p 'i-i " ^ s i +s i p ) p 'i + r i+i s i+i p 'i+i] 
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where s. - T / T. etc., are dimensionless quantities of order 1. 

i m i 

We see that T is a measure of the relaxation time. If t << t , the 
m me 

two terms on the left hand side of Eqs. (15) are vanishingly small compared to any 
term on the right hand side, since the magnitudes of the two terms inside the brackets 
on the left hand side are of order I . In this case, component 1 of the gas mixture 
in the i th interstice is "quasi -steady", which means that at any instant component 
1 in interstice i can be considered instantaneously steady and the derivative terms 
in its governing equation can be dropped. This argument applies to all other com- 
ponents of the mixture in every interstice, i .e., applies equally well to each 
equation of Eqs. (8) for cases j ^1 and Eqs. (9) individually. 

Asa final step towards completion of this argument, let 1 < m < L and 
1 < n < M be two integers such that T m ^ is the largest of all t ^ defined by 


,^= ( 4 V. /a. ) (tt m G)/8 k T.) 


1/2 


i i 


06 ) 
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Then, the relaxation time of the gas mixture in the entire insulation system T 
/ \ 

is bounded by L T ' ' , i .e., 
m 


T < LT W 
s ~ m 


07 ) 


Denoting the time period over which the properties of the gas mixture change 
appreciably by T , we can drop all the derivative terms in Eqs. (8) and (9) 
simultaneously if 

T or L T << t . (18) 

s m t 

Thus, this system of differential equations reduces to a set of algebraic 

.equations. . ^ 

_2 

Now, for ordinary engineering applications, t is of the order 10 ~ 

-1 s 
10 sec., while is measured by hours. Therefore, condition (18) is usually 

fulfilled, and Eqs. (8) and (9) reduce to 



i = 1, 2, , L " I- (20) 
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As Eqs. (19) and (20) are algebraic equations, there is no need for the pre- 
scription of initial conditions. 

To obtain the solution for a particular problem, we first choose a sequence 
of the time (t^, , t^). The solution at any chosen time, say t^, is ob- 

tained by solving Eqs. (19) and (20) with the chamber conditions and fuel tank 
temperature at t^ (see Eqs. (10) ) and (T, t^) (see Eq. (12) ) which are 
pertinent to that problem. After we finish with the time sequence, we have the 
pressures and temperatures in all the interstices of the insulation system as functions 
of the time. If the process is isothermal, one only has to solve Eqs. (19) with the 
relevant chamber conditions and outgassing rates. 
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4. Numerical Results and Comparison with Experimental Da'a 


For the purpose of comparison with experiment, we choose examples,the 
experimental data of which are available. 

A circular disk of 1 " thick and 6" in diameter is composed of a variety of 
numbers of insulation sheets of crinkled single-aluminized mylar. The edge of the 
disk is sealed with a solid insulation sheet (i.e., without perforation). The last 
layer (the back sheet) is also solid, but the rest of the insulation layers are per- 
forated with a perforation coefficient equal to 0.0138. The perforation holes are 
0.09375" in diameter. The distance between two adjacent holes is 0.707" center 
to center. 

The disk is placed in a large vacuum chamber. The chamber and hence 

the multilayer system (i.e., the disk) are maintained at room temperature all the 

time, so it is an isothermal process. At the beginning of the experiment, the 

vacuum chamber and the multilayer system are filled with nitrogen (purge gas), and 

the pressure inside each interstitial space between two successive insulation layers 

and the pressure in the chamber are atmospheric. The chamber is then evacuated. 

The pressure in the chamber and the pressure in the last interstice of the insulation 

system are recorded as time proceeds. 

/• \ 

The outgassing rates, R^'as functions of the time, of the insulation sheets 

were obtained before hand by separate experiments. The outgas components are 

mainly water vapor, with small amounts of nitrogen and carbon dioxide. For sim- 

( 2 ) 

plicity, we assume only one kind of outgas (water vapor) present, i.e., only R v 

( 2 ) 

has non-zero values. The outgassing rate R v ' (t) of crinkled single-aluminized 
mylar is shown in Fig . 1 . 

( 2 ) 

In our analysis, we use a set of R 7 (t) values at different times successively. 

(2) 

To obtain the solution at any chosen time, say t.,, we use the values of R v '(t j) 


0 ) 


and the chamber pressure p 7 (t .) in Eqs.( 19). Thus an appropriately chosen set 

(2) (1) cl 

of R' 7 (t) and p ' 7 (t) will result in a solution of the P. (t), the pressure history of 

C 1 2 

the gas mixture inside the multilayer insulation system. However, the R (t) values 

are available only up to t = 48 hours, thus our computations have to stop there. 
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Theoretical computation includes the pressures P.(t) inside each interstice, 
while only the pressure in the last interstice P^(t) , i .e., the backside pressure, is 
experimentally measured. Therefore, only the comparison for P^ (t) can be made. 
Comparisons of the theoretical P^ (t) and the experimental data are presented in 
Figures 2 to 6. 
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5. Conclusions 


First of all, it should be noted that the comparison between the computed 

and measured values of the p^'s presented in Figs. 2 to 6 cannot be taken seriously, 

since it is known that the experiment for the determination of the outgassing rate 

( 2 ) 

was not appropriately performed, thus the values of R v 1 (t) presented in Fig. 1, 

$ 

on which our calculation was based, is not reliable « Nevertheless, Figs. 2 to 6 
serve to show a qualitative comparison and to indicate some possible experimental 
errors . 

It can be seen from the figures that ar early limes the computed p^ is 
much higher than the experimental values. However, the computed and the mea- 
sured values are coming closer and closer toward each other when time proceeds, 
and good agreement is established after reasonably large time. Since the effect 
of preconditioning on the outgassing rate is likely to be a weak function of time 
when t is large, one probable and important course of experimental discrepancy 

would arise from different preconditionings of samples, i .e., the samples for the 

( 2 ) 

measurement of p^ (t)and the sample for the determination of R v '(t) were dif- 
ferently pre-treated. Another source of error would be the contamination of 
equipments inside the vacuum chamber. 

It is believed that agreement between the theory and experiment can be 
achieved if, with special attention to preconditioning, a set of consistent ex- 
periments is performed. 


At the time of the writing of this Report, the Marshall Space Flight Center is 
planning to reconduct the outgassing rate experiment. 
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FIG. 3. Comparison of Computed and Measured Backside Pressures, 20-Layer 
Sample. 
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n torrs 


c. 



FIG. 4. Comparison of Computed and Measured Backside Pressures, 35-Layer 
' Sample. . 
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